Duality of liquids 
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Liquids flow, making them remarkably distinct from solids and close to gases. At the same time, 
interactions in liquids are strong as in solids. The combination of these two properties is believed 
to be the ultimate obstacle to constructing a general theory of liquids. Here, we adopt a new 
approach to liquids: instead of focusing on the problem of strong interactions, we zero in on the 
relative contributions of vibrational and diffusional motion in liquids. We subsequently show that 
from the point of view of thermodynamics, liquid energy and specific heat are given, to a very 
good approximation, by their vibrational contributions as in solids, for relaxation times spanning 15 
orders of magnitude. We therefore find that liquids show an interesting duality not hitherto known: 
they are close to solids from the thermodynamical point of view and to gases from the point of view 
of flow. We discuss the experimental implications of this approach. 



The development of basic theories of solids involved a 
number of important discoveries that date back to over 
100 years ago. This was preceded by the development 
of the theories of gases. These developments form the 
basis of current understanding of most essential physical 
properties of these two basic states of matter. On the 
other hand, the third state of matter, the liquid state, 
remains poorly understood in comparison. A testament 
to this comes from the fact that even more recent text- 
books dedicated to liquids do not discuss most basic liq- 
uid properties such as specific heat [H. 

Liquids flow, and share this property with gases. At 
the same time, interactions in a liquid are strong, and are 
similar to those in solids. This presents a fundamental 
difficulty in calculating liquid energy in general form. In- 
deed, strong interactions, combined with system-specific 
form of interactions, imply that the energy is strongly 
system-dependent, precluding the calculation of energy 
in general form, contrary to solids or gases Q. 

Strong interactions are successfully treated in solids 
in the phonon approach, but this approach is thought 
to be inapplicable to liquids where atomic displacements 
are large. Stated differently, the "small parameter" in 
the theory of solids are atomic displacements, and a har- 
monic contribution to the energy, the phonon energy, is 
often a good approximation. The small parameter in 
gases are weak interatomic interactions. On the other 
hand, liquids have none of these because interactions are 
strong and displacements are large. The absence of a 
small parameter was, in Landau view, the inherent prop- 
erty of liquids, and was proposed as an explanation for 
why a theory of liquids can not be constructed at the 
same fundamental level as in solids or gases [1]. 

Here, we propose that stating a different reason for 
the difficulty of understanding liquids and exploring it 
provides an important way in to solving the problem. 
This reason is not unrelated to the one stated by Landau 
above, but operates at a different level. Instead of start- 
ing at the most fundamental level of strong interactions 
at the atomistic scale, we focus on the atomic trajectories 



that result from these interactions. 

We have recently proposed Q that from the point of 
atomic dynamics, solids and gases are pure states of mat- 
ter in the sense that dynamics in solids are purely os- 
cillatory and dynamics in gases are purely ballistic and 
coUisional. The dynamics of a liquid, on the other hand, 
is not pure but mixed: it involves both oscillations and 
ballistic motions [3] , and the relative contributions of the 
two types of motion change in response to external pa- 
rameters, temperature and pressure. Physically, the dif- 
ferent behaviors arise because in solids the kinetic energy 
of particles, K, is much smaller than the energy barri- 
ers between various potential minima, U: K <^ U. In 
gases, it is the other way around: K :s> U. In liquids, 
the mixed nature of dynamics originates because none of 
these strong inequalities apply. 

Notably, the solid-like oscillatory component of liquid 
dynamics originates from large energy barriers preventing 
local diffusion processes. Large energy barriers are set 
by strong interactions - this is how interaction strength, 
emphasized by Landau, enters our approach. 

In our approach to liquid thermodynamics, we there- 
fore focus on atomic trajectories and relative weights of 
solid-like oscillatory and gas-like diffusional motions. In 
this approach, we demonstrate close similarities of liquid 
and solid thermodynamics, and consequently propose a 
new property of liquids, their duality. Our next step is 
to calculate the liquid energy by evaluating relative con- 
tributions of oscillatory and ballistic diffusional motions. 

The evaluation can be done using the concept of relax- 
ation time proposed by Frenkel Atoms or molecules 
in a liquid are not fixed, but rearrange in space due to 
thermal fluctuations, giving liquid flow. Each flow event 
is a jump of an atom from its surrounding cage, accompa- 
nied by large-scale rearrangement of the cage atoms. We 
call this process a local relaxation event (LRE). A LRE 
lasts on the order of Debye vibration period td = 0.1 ps. 
Frenkel introduced liquid relaxation time r as the time 
between LREs at one point in space in a liquid r is 
directly related to liquid viscosity rj via the relationship 
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V = GoqT, where Goo is the instantaneous shear modu- 
kis. At high temperature, r reaches its hmiting value, 
T « td, when time between the jumps is equal to the 
shortest vibrational period. On temperature decrease, r 
increases until it reaches the value at which the flow stops 
at the experimental time scale. This corresponds to the 
liquid-glass transition and r w 10"^ s. 

Lets consider the regime where LREs take place rarely 
compared to the short period of vibrations: 
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time T, all LREs in the system relax. Therefore, we need 
to calculate ii^dif that is averaged over time r. 

Let us divide time r into m time periods of duration 
td each, so that m 
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Then, i?dif , averaged over time 
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where E^-^^ are instantaneous values of i?dif featured in 
Eq. ©. g is 



There are two contributions to total liquid energy, Ei: 
El = Evih + Edif, where E^ih and i?dif are the energies 
of vibrations and diffusion, respectively. -Evib m a liq- 
uid is given by one longitudinal mode and two transverse 
modes with frequency cj > ^ (a liquid docs not support 
transverse modes at lower frequency, but yields to shear 
stress instead) When the regime ([Ij is operative, 
-Evib to a very good approximation is given by 3A^ modes, 
Es — 3NT, as in a classical harmonic solid (here, N is the 
number of atoms and fee = 1). Indeed, oj > - can be re- 
stated as a; > ujf, where is Frenkel frequency. When 
T ^ Td, wp <C wd, where = is Debye frequency. 
This implies that most of the transverse modes existing 
in a solid are supported in the liquid when r ^ td. Fur- 
thermore, it is only the high-frequency modes that make 
a significant contribution to the liquid vibrational energy, 
because the density of states of quasi-harmonic modes is 
approximately proportional to o;^. 

Edii is due to atoms involved in LREs, and includes the 
kinetic energy of jumping atoms and the energy of their 
interaction during the jump events. The jump probabil- 
ity for one LRE, p, is the ratio between the time spent 
diffusing and vibrating: p = In statistical equilib- 
rium, p is equal to the ratio of diffusing atoms, iVdif, 
and the total number of atoms, TV. Then, at any given 
moment of time: 
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If -Edif is the energy associated with diffusing LREs, 
i?dif oc A'^dif- Together with Etot oc N, Eq. ^ gives 
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Eq. ([3]) implies that under condition ((T)) , the contribu- 
tion of _Edif to the total energy at any moment of time 
is negligible. This time is set by the smallest time scale 
of the system, td, and corresponds to the instantaneous 
value of i?dif • During time time, the system in not equi- 
librium. The equilibrium state is reached when the ob- 
servation time exceeds system relaxation time, r. After 
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Each of the terms in Eq. ([5]) is equal to ac- 
cording to Eq. ([3]) . There are m terms in the sum in Eq. 
(H. Therefore, 
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We therefore find that under the condition ([T]), the 
ratio of the average energy of diffusion motion to the 
total energy is negligibly small, as in the instantaneous 
case, Eq. ([3|). Consequently, the energy of the hquid 
under the condition ([T]) is, to a very good approximation, 
given by the remaining vibrational part. As discussed 
above, the vibrational part of the liquid energy under 
the condition ([IJ is, to a very good approximation, given 
by the vibrational energy of all 3iV modes as in a solid. 
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This, in turn, implies that the liquid constant-volume 
it, c 

regime Eq. ([1]): 



specific heat, Cyj — is entirely vibrational in the 



Cv,l 



(8) 



The primary evidence supporting Eq. ([5]) and our the- 
ory comes from experimental Cy^i. Experimentally, Cy^i of 
many systems including metallic, noble, molecular and 
network liquids is very close to 3, the Dulong-Petit value 
of a solid, when Eq. ^ holds 0-1^. For low- viscous liq- 
uids this corresponds to temperatures around the melting 
point. 

We note that on temperature increase when condition 
T starts to approach td and Eq. ([T]) no longer applies, 
experimental Cy^i starts to decrease from about 3 at the 
melting point to 2 at high temperature 0-01, an effect 
that we recently attributed to the progressive loss of 
transverse modes with frequency w > ^ (8|. 
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As discussed above, the vibrational energy and heat 
capacity of the liquid under the condition ^ are given by 
3A'' vibrational modes. In the harmonic case, this gives 
the energy and heat capacity of a liquid as 3NT and 3A^, 
respectively, i.e. the same as in solids. Then, Eqs. (|7l8p 
imply the equality of liquid and solid vibrational energies 
and specific heats in the harmonic case: 



Ej^ = E^ 



(9) 



(10) 



where s corresponds to the solid and h to the harmonic 
case. 

In the anharmonic case, Eqs. ([7]) and Eqs. ^ still 
hold, but Eqs. © and Eqs. pUj) hold only approximately 
because anharmonicity affects the vibrational energy and 
Cv of a liquid and a solid in a different way. In particular, 
Cv is modified by the intrinsic anharmonicity related to 
softening of vibrational frequencies at constant volume, 
and can be approximately written as c„ = 3(1 + aT), 
where a is the coefficient of thermal expansion a is 
generally larger in liquids compared to solids, resulting in 
larger Cv,i compared to Cv,s, albeit the term aT is usually 
small compared to 1. 

Interestingly, although Eq. ([6]), combined with Eq. 
P^, implies that c^ of a liquid is entirely vibrational 
as in a solid, this does not apply to entropy (here and 
below we imply the equilibrium state, the condition for 
which is t ^ r, where t is observation time): the dif- 
fusional component to entropy is substantial, and can 
not be neglected. Indeed, if Zvib and Z^^i are the con- 
tributions to the partition sum from vibrations and dif- 
fusion, respectively, the total partition sum of the liq- 
uid is Z = Zvib • Z^if. Then, the liquid energy ]s E — 
(ln(Zvib ■ Zdif)) = T^^lnZvib + T^^ln^dif = 
ii^vib + E^if (here and below, the derivatives are taken 

at constant volume). Next, ^ 1 from Eq. ^ also 

implies 1, where, for brevity, we dropped the sub- 

script referring to the average. Therefore, the smallness 
of diffusional energy, <^ 1, gives 
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dT 

The liquid entropy, S — (Tln(Zvib • ^dif))j is: 



S = T-^ In Zvib + In Zvib + T-^ In Zdif -t- In Zdif (12) 

The condition ()11|) implies that the third term in Eq. 
p2p is much smaller than the first one, and can be ne- 
glected, giving 
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Eq. ([T3| implies that the smallness of i?dif, expressed 
by Eq. (fTTj) . does not lead to the disappearance of all 
entropy terms that depend on diffusion because the term 
In Zdif remains. This term is responsible for the excess 
entropy of liquid over the solid. On the other hand, the 
smallness of E^a does lead to the disappearance of terms 
depending on Z^a in the specific heat. Indeed, cj, = T^y 
(here, S refers to entropy per atom or molecule), and 
from Eq. ()13|) . we find: 
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Using Eq. pT|) once again, we observe that the third 
term in Eq. p4p is small compared to the second term, 
and can be neglected, giving 
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dT V dT 
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As a result, c„ does not depend on Zdif, and is given 
by the vibrational terms that depend on Zvib only. As 
expected, Eq. P3)) is consistent with Eq. ([5]). 

Physically, the inequality of liquid and solid entropies. 
Si ^ Ss (despite Eq. ([TT]) ). is related to the fact that the 
entropy measures the total phase space available to the 
system, which is larger in the liquid due to the diffusional 
component present in Eq. (1131) . This diffusional compo- 
nent, although large, is slowly varying with temperature 
according to Eq. (jlip . resulting in small contribution to 
Cy (see Eqs. and ((T5|) ) and giving Cy,i = c^^s- As 

far as the energy is concerned, it reflects the instanta- 
neous (or averaged over r) state of the system that does 
not rely on exploring the phase space. Consequently, 
El = Eg, yielding Eq. ([TT|) and the smallness of diffu- 
sional contribution to Cy. 

We conclude that two important properties of a liquid, 
energy and specific heat, are essentially vibrational, as 
they are in a solid, provided ^ -C 1. For practical pur- 
poses, this takes place for r > IOtd. Perhaps not widely 
recognized, the condition t w IOtd holds even for low- 
viscous liquids such a liquid monatomic metals (Hg, Na, 
Rb and so on) and noble liquids such as Ar near their 
melting points 0, S], let alone for more viscous liquids 
such as room-temperature olive or motor oil. 

Notably, the condition t > IOtd corresponds to almost 
the entire range of t at which liquids exist. This fact 
was not fully appreciated in earlier theoretical work on 
liquids. Indeed, on lowering the temperature, r increases 
from its smallest limiting value of r = td ~ 0.1 ps to 
T « lO'^ s where, by definition, a liquid forms a glass at 
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the glass transition temperature. Here, r changes by 16 
orders of magnitude. Consequently, the condition ^ ^ 
1, Eq. or T > IOtd, applies in the range 10^ - 10"^^ 
s, spanning 15 orders of magnitude of r. 

In addition to explaining the experimentally observed 
Cv , our approach provides general insights into the prob- 
lem of liquid-glass transition, and implies that the theory 
of solids is a good starting point for discussing important 
aspects of liquid thermodynamics. This is a useful result 
in the area where the calculation of thermodynamic prop- 
erties of viscous liquids approaching glass transition has 
been a long-standing problem [lo| . 

We can also revisit the long-standing and puzzling 
question of why Phillips constraint theory of glasses |llj] 
works. Proposed over 30 years ago, the theory is based on 
the idea that a short-ranged interatomic bond in a glass 
can be viewed as a rigid mechanical constraint. Equating 
the number of constraints and degrees of freedom gives 
the average coordination number (r) —2A at which, the 
theory proposes, the glass-forming ability is optimized. 
The constraint theory of glasses has since enjoyed success 
in explaining other properties of glass transition [l^ . 

The constraint theory derives its predictions from 
counting the bonds (constraints) and (r) in the solid 
glass. Importantly, the theory subsequently uses (r) to 
predict the properties measured across the liquid-glass 
transition and in the liquid state. This is truly surpris- 
ing: indeed, in the liquid state above Tg, bonds are not 
intact because atoms rearrange on the experimental time 
scale, and therefore can no longer be viewed as rigid con- 
straints. This poses an intriguing question of why the 
constraint theory works. This question remained an open 
and long-standing puzzle in the area. 

Our approach to liquids readily answers this question 
on general grounds: the measured properties operate in 
the regime where Eq. ([T]) applies. Indeed, glass-forming 
melts such as silicates, chalcogenides and others are mea- 
sured at temperatures where they are viscous enough to 
satisfy Eq. ^ [l2|. Consequently, our theory predicts 
that if a measured property is related to system energy or 
specific heat, this property in the liquid state is, to a very 
good approximation, equal to that in the solid glass, and 



can therefore be predicted from solid-state properties. 

That important elements of liquid thermodynamics 
can be understood on the basis of thermodynamics of 
solids is a result hitherto not anticipated. Indeed, start- 
ing from earlier proposals [2], existing theories approach 
liquids as interacting gases and consequently attempt to 
calculate liquid energy as an integral of the product of in- 
teratomic interactions and correlation functions 1] . The 
interactions and correlations are often complex and are 
not generally known, except in the simplest liquids such 
as Ar. Our result circumvents this problem because the 
thermodynamics of solids is mostly understood. 

In summary, we observe that liquids flow, making them 
remarkably distinct from solids and close to gases. How- 
ever, from the point of view of thermodynamics, liquid 
energy and specific heat become, to a very good approxi- 
mation, equal to those in solids for relaxation times span- 
ning 15 orders of magnitude. In this sense, liquids show 
an interesting duality of physical properties. 
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